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We demonstrate a method to extend the range of pulsed laser spin noise measurements to long spin lifetimes. 
We use an analog time-domain detection scheme with a bandwidth limited only by laser pulse duration. Our 
model uses statistics and Bloch-Torrey equations to extract the Lande g-factor, Faraday cross-section ap, 
and spin lifetime rs, while accounting for finite detector response. Varying the magnetic field with a fixed 
probe-probe delay yields ts when it is longer than the laser repetition period. Varying the probe-probe delay 
with a fixed field produces a time-domain measurement of the correlation function. 


Noise is an expression of a system’s fundamental be¬ 
havior. The fluctuation-dissipation theorem^ii^ when 
combined with the Wiener-Khinchin theorem,^ implies 
that a system’s dynamic response can be extracted from 
correlation measurements. Spin noise was theoretically 
postulated in 1946^ and then demonstrated on atomic 
gases,— nuclear spins^ and with scanning tunneling 
microscopy— More recently, spin noise spectroscopy on 
semiconductor o^°d^ has revealed its potential for per¬ 
forming contactless and non-perturbative measurements 
that can also more directly access the Faraday cross- 
section. 12, Experiments that utilize spin noise as a level 
of contrast have performed non-destructive spatial reso¬ 
lution of a semiconductor’s dopant densityi^ and homo¬ 
geneous absorption linewidth measurements of quantum 
dot ensembles—^ 

Conventional spin noise spectroscopy utilizes a contin¬ 
uous wave laser transmitted through a sample to measure 
the Faraday rotation produced by stochastic spin fluctu¬ 
ations. A magnetic field is applied orthogonal to the 
optical path to induce precession and measure the Lande 
g-factor. A power spectrum is produced either by using 
electronic sampling and Fourier transforming the data 
or using a spectrum analyzer— Sidi Isolation of 
the spin noise contribution from a typically much larger 
background is accomplished by taking the difference of 
the finite power spectra for two system states. In these 
measurements, the bandwidth is limited by the sampling 
rate. 

In order to access high-frequency dynamics, 
theoreticali^ and experimentali2ri2i efforts have in¬ 
cluded the use of ultrafast pulsed lasers. Ultrafast 
spin noise spectroscopy allows for a direct time-domain 
measurement of the spin correlation function, with a 
bandwidth determined by the laser pulse duration / ^ i 22 
However, the temporal range is limited by the optical 
delay time between probe pulses and, ultimately, the 
laser repetition period tr-ep, which makes it difficult to 
accurately determine the spin lifetime ts when ts ^ Uep- 
In this Letter, we demonstrate a Resonant Spin Noise 
(RSN) technique that enables accurate measurements 
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of the spin lifetime when ts ^ trep and clarify how to 
theoretically model the ultrafast spin noise signal and 
extract the Faraday cross-section. We compare our data 
and model with the model presented in Ref. 21 which 
predicts a Lorentzian peak at zero magnetic field that 
we do not observe. 

We perform ultrafast spin noise measurements using 
analog electronics. We process our signal with an ana¬ 
log root-mean-square (RMS) circuit and then isolate the 
spin noise component with lock-in amplification. We de¬ 
velop a model to explain the analog electronic calculation 
based on statistics and simplified Bloch-Torrey equations 
which are derivable from the semi-classical Boltzmann 
equation— i22i21 By including the response time of our cir¬ 
cuitry, we are able to quantify our signal amplitude and 
extract the Faraday cross-section— 2 We show that when 
the spin lifetime ts is short compared to the laser repe¬ 
tition period trep, measurements as a function of relative 
probe delay At accurately capture the spin dynamics. 
When Ts ^ trep, we show that measurements as a func¬ 
tion of applied magnetic field with fixed At accurately 
determine ts- 

A system’s dynamics can be obtained from noise 
through the correlation of two measurements 9i and 9j. 
In practice, the correlation function is obtained from 
measurements of the expectation value of the product 
of 9i and 9j: 

{9,9,) = (9,) (9,) + Var {9) Corr {9„ 9,) (1) 

Here, we separate the correlation into the product of its 
magnitude, Var [9), the variance of measurements 9i, 
and a normalized correlation function Corr (9i,9j). For 
a sample with no net spin polarization, (0/ = {9j) = 0. 

To acquire {9i9j), we can sum 9i and 9,, square the 
sum, and then perform an average of multiple squared 
sums. We accomplish these operations using a series of 
analog electronic components including a low pass fil¬ 
ter and the RMS channel of a telecommunications chip 
(Analog Devices ADL5511-EVALZ). Due to the analog 
processing, the first sum turns out to be a weighted aver¬ 
age (WA) of many 9i, as explained in the Supplementary 
Material.^! To isolate the spin contribution to the signal 
9 from the background we use lock-in detection, with 
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the result 

RMS Signal = ((0 + " {^WA?'^ 

— {^wa) ^ ( 2 ) 


In the second line, we have made the assumption that 
0 ^ ^ which, for spin noise measurements of n- 
GaAs, can be accomplished through the use of sufficient 
probe power, although too much power will perturb the 
samplej^ (^^a) is defined as: 
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The leading integral averages the voltage decay between 
the most recent probe pulse and the next probe pulse 
where tj^c is the electronic response time. The discrete 
sums provide a weighted average over all previous pulse 
residues up to, and including, the most recent pulse. 

Our model^^ for carrier spin dependent behavior in n- 
GaAs is 


(OiOj) = Var (0) e — tj\) (4) 

with Var {9) = apnd/16A where n is the sample carrier 
density, d is the sample thickness, and A is the probe 
cross-section at the sample, ctf is the Faraday cross- 
section where 9 = apnpd with ns = nf — n^ = Pn 
being the sample’s spin polarized carrier density with P 
the level of polarization.^^ 12 = gBgB/h is the Larmor 
precession frequency where is the Bohr magneton, g 
is the Lande g-factor, h is the reduced Planck’s constant, 
and B is the applied magnetic field strength. 

We interpret the Bloch-Torrey equations as deriving 
the ensemble correlation function, a normalized descrip¬ 
tion of the dynamics. Our function Var (9) contains all 
information about the amplitude and is calculated by 
taking an RMS over all initial ensemble states. As the 
Bloch-Torrey equations are semi-classical and vectorial, 
we represent the two-state system for electrons with a 
binomial distribution. We also average over all possible 
orientations on the unit sphere with the assumption that 
there is no preferential axis of quantization. 

Our experiment utilizes a 76 MHz repetition rate 
{trep — 13.6 ns) Ti:Sapphire laser with 3 ps pulse dura¬ 
tion. We pass the output through a beam splitter (BS) to 
create fixed and mechanically variable length (VL) paths 
with separate pulse trains. We recombine the paths at 
a second beam splitter prior to entering a fiber optic ca¬ 
ble (FOG). After exiting the FOG and passing through 
a linear polarizer, the two pulse trains will be identical 
in polarization, intensity, and spatial profile prior to in¬ 
teraction with the sample where they undergo Faraday 
rotation proportional to the stochastic spin polarization. 
Each pulse train has an average power of 10 mW and is 
focused to an area of 500 gvc? at the sample. 

We use an electro-optic modulator (EOM) to turn on 
and off our ability to detect the Faraday rotation of pulses 


transmitted through the sample, as described in Ref. 21. 
This is achieved by either allowing unaltered linear polar¬ 
ization to proceed (EOM is at zero retardance) or chang¬ 
ing the linear polarization to circular polarization (EOM 
is at quarter-wave retardance) with a 157 Hz square-wave 
switching frequency. After the EOM, the light passes 
through a Wollaston prism (WP) which splits the pulses 
into orthogonal linear components. When the EOM is 
at zero retardance, a Faraday rotation will be detected 
as an imbalance at the photodiode bridge (PDB). At 
quarter-wave retardance, no imbalance will be detected. 
After the PDB, the photogenerated voltage is amplified, 
and then processed on the RMS circuit of the ADL5511- 
EVALZ. The limited speed of our PDB, combined with 
parasitic capacitances in our wiring, provide an effec¬ 
tive low-pass filter to generate a weighted average with 
tpc 29 ns. Final detection of our modulated signal 
occurs at a lock-in amplifier (LIA). 

We tested our measurement scheme by performing 
spin noise measurements on a 500 /rm thick chip of 
bulk Si-doped n-type (001) GaAs. The following are 
the manufacturer specifications: carrier concentration of 
(4.3 —6.2) X 10^® cm“^, mobility of (3450-3880) cm^/V-s, 
and resistivity of (2.8—3.9) x 10“^ Ohm-cm. The electron 
carrier concentration of ~ 5x 10^® cm“^ is near the metal- 
to-insulator transition of GaAs where long spin lifetimes 
are expectedi^ 

Figures [Hand [5] show Resonant Spin Noise (RSN), in 
which the delay between two pulse trains remains fixed 
while the magnetic field is swept. The resultant signal is 
reminiscent of the pump-probe method of Resonant Spin 
Amplification (RSA)^ but arises from a different mea¬ 
surement principle. RSA is the result of multiple induced 
spin populations with integer phase relationships under¬ 
going precession that are all probed at the same point in 
time. RSN is the precession of a single stochastic spin 
population, probed at a series of times that undergo sub¬ 
sequent processing to extract correlations. Due to its 
minimally perturbative capabilities, RSN, and spin noise 
more generally, should provide simpler analysis of ts 
compared to values extracted with RSA if the optically- 
pumped spin polarization is affected by subsequent pump 
pulses, as observed in Ref. 27. 

We compare data collected at 10 K for a single pulse 
train and two pulse trains with different delays. Using 
a relative delay (At ^ 6.6 ns = Uep/2) that is half of 
the laser repetition period leads to a signal (Fig. |T1 top) 
that is similar to a single pulse train (Fig. [H bottom). 
This is due to the equal time spacing of the probe pulses. 
Equations [H |3l and [4] are used to model the behavior of 
the single pulse train data, shown as a solid gray line, 
with <5 = trep and ti — tj = (i — j)trep- Letting t^ep —>■ 
tre-p!^ allows for fitting of the two pulse train case. 

For probe pulses that are not equally spaced, we must 
modify Eq. [31 Equation |S] is a weighted average of the 
two possible timing sequences that stem from the first-in, 
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FIG. 1. Resonant Spin Noise (RSN) data taken at 10 K with fits shown as solid gray lines. Single pulse train data is fit using 
Eqs. [HO and [4] while two pulse train data is fit by Eqs. [4] and (5) Data are scaled and offset for clarity. 



FIG. 2. Resonant Spin Noise (RSN) data taken at various temperatures with fits using Eqs. [4] ^nd [5] shown as solid gray lines. 
I At I = 2 ns for all data. Data are scaled and offset for clarity. 


first-out behavior of our electronic processing. 

{^Wa) = \/ fa,P + - \/ fp,a (5a) 

Lrep ^rep 

(5b) 

with tf G {0,a,trep,trep + ...}, tt = |At|, and /3 = 

trep — I At I where |At| < trep is the relative spacing be¬ 
tween the two pulse trains controlled by the variable 
length optical delay line. Equations 0] and [S] are used 
to fit all two pulse train data, including Time-Resolved 
Spin Noise (TRSN) shown in Fig. [3] We calibrated our 
system using |At| = 2 ns RSN data, combined with g 
and Ts values extracted from 10 K RSA measurements, 
yielding tuc ^ 29 ns. It should be noted that this leaves 
only three parameters for the rest of the fits: g, ts, and 
(jp- 

Figure 0] shows the temperature dependent evolution 
of RSN data with |At| = 2 ns up to 75 K. As temper¬ 
ature increases, ts decreases as well as the contribution 
of terms with larger temporal spacing in Eq. [5] By 75 
K, rs is small enough that the error in fitting RSN is ap¬ 


preciable. However, the 75 K data is a useful qualitative 
test of our model. We do not observe a Lorentzian peak 
at zero magnetic field as expected from the antiderivative 
of Eq. I in Ref. 21^ 

To obtain a more accurate measurement of the spin 
lifetime for shorter ts, we can instead perform TRSN 
where we scan the time delay At between the pulses at 
a fixed magnetic field. In Time-Resolved Faraday Rota¬ 
tion (TRFR)^ a pump pulse orients a single population 
of spins and a delayed probe measures their evolved ori¬ 
entation. In TRSN, a single population of stochastically 
aligned spins is probed by an initial pulse and then com¬ 
pared to a delayed pulse through Eq. 01 capturing the 
systems evolution. 

Figure 0] shows TRSN at 100, 150 and 200 K with 
an applied magnetic field of 300 mT. Equations 0] and 
|S] are used for fitting and shown as solid black lines. If 
|At|,Ts <C tsc,trep, we Can simplify Eqs. 0] and |S] to 
obtain Eq. 01 

^ \J^ + e-l^*l/^scos (H|Af|) (6) 

Figure m shows the temperature dependence of the fit 
values g and ts determined by RSA and TRFR and our 
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FIG. 3. Time-Resolved Spin Noise (TRSN) data taken at 
various temperatures with fits using Eqs. |4] and [5] shown as 
solid black lines. The applied magnetic field strength was 300 
mT for all data shown. Data are scaled and offset for clarity. 



Temperature (K) 

FIG. 4. Temperature dependence of fit values extracted from 
Resonant Spin Noise (RSN) and Time Resolved Spin Noise 
(TRSN) data shown in Figs. [3]and|4]compared with fit values 
extracted from the conventional techniques of Resonant Spin 
Amplification (RSA) and Time Resolved Faraday Rotation 
(TRFR). The blue dashed line is a power law fit, ^ ^ 

for the D’yakonov-Perel dephasing mechanism. 


spin noise methods of RSN and TRSN. The agreement is 
quite good, with the exception of the RSN value for ts 
at 75 K, and is consistent with previous reports i^^'^^i^° 
We have included a blue dashed line in Fig. 0] showing 
the expected D’yakanov-Perel temperature dependence 
for Ts which, for T > 50 K, follows the power law rg ~ 

jn-5/2^ 

Figure El shows the temperature dependence of 
\Var (0)]^^^ extracted from RSN and TRSN data. Since 
our sample is near the metal-to-insulator transition and 
the energy level of Si donors in GaAs is ~ 6 meV below 
the conduction band(^ the majority of Si donors should 
be ionized above 100 K. We observe nearly constant be¬ 
havior for [Rar(0)]^^^ above 100 K and attribute this 
to a constant carrier density, which results in a constant 



FIG. 5. Temperature dependence of Var{6)^^^ extracted 
from RSN and TRSN data. The red dashed line is the aver¬ 
age of Var for 100 K through 200 K used to extract a 

value for ap. 

value for ap. We use the average value of [Var {9)]^^'^ 
above 100 K, along with n = 5 x 10^® cm“®, d = 500 
/im, and A = 500 /rm^ to extract ap yielding the value 
ap ~ lx 10“^^ rad-cm^ in agreement with Ref. 12. More¬ 
over, our measurement does not require any assumptions 
regarding optical pumping efficiency. 

Spin noise techniques can also offer advantages when 
attempting to study material systems with inefficient op¬ 
tical orientation. The signal strength of conventional 
pump-probe techniques relies upon spin dependent op¬ 
tical selection rules, as well as sample dimensions and 
probe spot size. We can determine the approximate 
threshold where spin noise provides a stronger signal 
through the ratio of \Var (0)] ' to the maximum pump- 
probe Faraday rotation Opp = apPnd: 

dpp APy/nAA 

The dimensionless maximum polarization value P is de¬ 
termined by selection rule efficiency and experiment con¬ 
ditions. If P < l/A\/ndA, then spin noise should be a 
more advantageous experiment. We can use our result 
to evaluate the interesting case of probing a single spin 
where ndA = 1. Spin noise should allow the observation 
of dynamics, without the need to filter out pump scatter, 
with signal that is on the order of maximum achievable 
pump-probe Faraday rotation. 

In conclusion, we have demonstrated Resonant Spin 
Noise (RSN) and Time Resolved Spin Noise (TRSN). 
Both methods utilize the same optical path and elec¬ 
tronics which, in principle, offers a combined bandwidth 
beyond 1 THz with commercially available femtosecond 
pulsed lasers. We model our analog processing and 
explored spin dynamics using statistics and simplified 
Bloch-Torrey equations. We interpret the Bloch-Torrey 
equations as deriving the correlation function while the 
amplitude of the signal is the root-mean-square of all pos¬ 
sible system states. From fitting, we extract values for 
the dephasing time, Lande g-factor, and Faraday cross- 
section that agree with conventional measurements and 
available literature. 
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FIG. S 1. Optical path described in the text. BS = beam 
splitter, EM = electromagnet, EOM = electro-optic modu¬ 
lator, FOG = fiber optic cable, LIA = lock-in amplifier, LP 
= linear polarizer, PDB = photodiode bridge, VL = variable 
length path, WP = Wollaston Prism. 


Appendix A: Derivation of {6i9j) for spin dependent 
behavior in n-GaAs 

We describe the measured spin physics of our sample 
using Eg [AH 

^-nxS + — = SoS{t) (Al) 

ot TS 

where S is the time-dependent spin density, So is the ini¬ 
tial spin density, $7 is the Larmor precession frequency 
as defined in the Letter, ts is the spin dephasing time, 
and 6 (t) is the Dirac-delta function. The difference be¬ 
tween the spin up and spin down components of the 
carrier density n along the optical path is defined as 
ns = n-j- — 714 , = |So| with n = + ni. Equation lAll is 

an isotropic version of the Bloch-Torrey equations that 
excludes drift and diffusion and is applicable to our n- 
GaAs sample and experiment conditions. We have in¬ 
corporated 6 (t) so that we may solve for the Green’s 
function. 

Using a sphere (see Eig. we define the spin projec¬ 
tion along the optical path to be in the x-direction with 
an applied magnetic field in the z-direction yielding: 

Sx = nsH [t — to) sinc/ig cos [U (t — to)] (A2) 

where (jjo and /3o are the initial orientation of the spin 
polarization. H {t — to) is the Heaviside step function 
with to being the initial time of generation for the spin 
population. 

We let t — to —>■ ti and Sx —t Sxi be a random time and 
associated x-projection of the spin polarization respec¬ 
tively. We also ignore the Heaviside function by assum¬ 
ing repeated sampling of a single spin population that 
was polarized at a random instant in time. We then av¬ 
erage over all configurations of ns, Po, and 4>o- As in Eq. 
[11 we obtain the correlation from the expectation of the 
product SxzSxj 


{SxiSxj) = (Sxi) (Sxj) + Var (Sx) Corr {Sxi, Sxj) (A3) 
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I (down) 


FIG. S 2. Sphere for the vectorial representation of the spin 
density vector. 


For a sample with zero net spin polarization, {S^i) = 
{Sxj) = 0, and {Sxi) (Sxj) = 0. We write the expectation 
value of the product of two ensemble spin states as: 

{Sxi^xj') — 

sin^ 4>o cos [flti + /3o] cos [fltj + /3o]^ 

(A4) 

Only the relative time, not the absolute time, between 
states is of importance so we let tj ti + |At|. We also 
let —>■ 0 as the absolute time should not affect the value 
of the correlation between measurements, yielding: 

(SxiSxj) = sin^ (j)o cos/3o cos [fl\At\ + /3o]^ 

(A5) 

We can relate the measured Faraday rotation to the 
ensemble spin polarization using 0 = upnsd. We can 
also simplify our expression by evaluating the average of 
(j)o G [0,7r) and /3o S [0,27r) which results in the factor 
1/4. Our result is 

{OiOj) = i {a-pd)'^ (^1) cos (0| At|) (A6) 

where we recognize Var{6i) = j{(jFd)^ (rig'j and 
Corr{9i,9j) = As (.Qg At|). 

The last portion of our derivation involves the evalua¬ 
tion of the radial magnitude of the spin polariza¬ 

tion. Using the fact that spins are a two state system, 
we can evaluate (ug) as the variance of a binomial dis¬ 
tribution. Therefore (n|) = N/AV'^ where N is the total 
number of probed spins and V = Ad is the probed volume 
equivalent to the probe cross-section multiplied by the 
thickness of the sample. If we recognize that n = N/Ad, 
we arrive at our hnal expression for the variance: 


Var {9^) 


ap nd 

16 T 


(A7) 


Appendix B: Discussion of weighted average for two pulse 
trains with variable delay 

We begin with a schematic representation of the elec¬ 
trical processing. On the left in Fig. 3S1 the yellow boxes 
highlight the two possible scenarios: summing over all 
previous pulses while averaging the decay over an inter¬ 
val of At (top, series in blue); summing over all previous 
pulses while averaging over the decay over an interval of 
trep — At (bottom, series in red). We treat the electrical 
pulses generated at the photodiode bridge as having an 
instantaneous rise time (3 ps Ti:Saph laser pulsewidth) 
and a finite fall time due to the limited response of the 
electronics {tpic 29 ns). 



t„p - At Weiglitecl Simi 



FIG. S 3. Schematic representation of electrical processing. 

Since our electronics are in series, the processing is 
first-in-first-out which we schematically follow on the 
right side of Fig. 321 First, a weighted sum over ei¬ 
ther the blue series or the red series is performed. Sec¬ 
ond, the value of the weighted sum is squared and then 
square-rooted. Third, an average over many weighted 
sums is performed leading to a steady state value which 
we extract using lock-in detection. 

We can write down the mathematics in a straightfor¬ 
ward fashion. For the squared pulse response weighted 
sum, we write 


/«./3 


1 



2 






(Bl) 

Equation IBlI is a modified version of Eq. 3 in the pa¬ 
per where the time-points are not evenly spaced. For 
/a./ 3 , /3 is the period for calculating a given weighted sum 
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FIG. S 4. Plot of combined Eqs. IB1IIB2I and IClI for the same 
parameters as in Fig. [2] in our paper. Equation ICll comes 
from taking the antiderivative of Eq. 1 in Ref. 21. As shown 
in our paper, we do not observe a pronounced peak at 0-field 
in our data set. 


while a is the first non-zero time-point in the series of 
sampled times G {0, a, t^ep, trep-i-a, }■ For our experi¬ 
ment, /? = At while a = t^ep — At with correspond¬ 
ing to the blue series and fp^a (the permutation of a and 
/3 values) to the red series shown in Fig. 331 
As shown on the bottom right hand side of Fig. 331 we 
only measure a single value as our electronics perform a 
square-root, followed by a weighted average, of fa,p and 
fp,a- We write this mathematically as: 

{^Wa) = ~. \/ fa,P + - \/ fp,a (B2) 

Lrep ^rep 


where we make use of the fact that a + {3 = trep- The 
weighting factor for each series is determined by the rel¬ 
ative interval lengths of the respective weighted sums. 


Appendix C: Comparison to Model Used by Ref. 21 

If we take the antiderivative of Eq. 1 in Ref. 21, we 
arrive at: 


(0 0 ) ^ _ 1 g-|At|/Ts 

X ((2 -I-T|n^) cos [n|At|] — rsfl sin [n|At|]) (Cl) 

This would predict a Lorentzian type scaling near zero 
magnetic field for a Resonant Spin Noise scan. If we use 
this model for Eqs. IBII and IB2[ we would expect the 
signal plotted in Eig. 331 

Reference 21 does not include a full derivation of their 
model. However, Eq. EH can be derived by calculating 
the integral 


dtfit)fit-\At\) (C2) 

/ (t) = cos [Ht] (C3) 

Equation EH is the autocorrelation of f (t), but the re¬ 
sult (Eq. ICll) does not appear to accurately describe the 
measured spin noise behavior. In our model (Eq. IA6I) . 
the spin ensemble correlation function is Eq. IC3l and not 
its autocorrelation. 








